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Abstract 

We explore the physics potential of the NuMI beamline with a detector located 10 km off-axis 
at a distant site (810 km). We study the sensitivity to sin 2 2#i3 and to the CP-violating parameter 
sin 5 as well as the determination of the neutrino mass hierarchy by exploiting the — ► v e and 
Vfj, — ► v e appearance channels. The results are illustrated for three different experimental setups to 
quantify the benefits of increased detector sizes, proton luminosities and v e detection efficiencies. 
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Neutrino oscillations have been observed and robustly established by the data from so- 
lar □ D atmospheric jjj, reactor jj| and long-baseline neutrino experiments j^. These 
results indicate the existence of non-zero neutrino masses and mixings. The new parameters 
can be accommodated via the three neutrino PMNS mixing matrix 1 , the leptonic analogue 
to the CKM matrix in the quark sector. Neutrino oscillations within this scenario are de- 
scribed by six parameters: two mass squared differences 2 (Am^ and Am^), three Euler 
angles (#12, #23 and #13) and one Dirac CP phase S. The standard way to connect the solar, 
atmospheric, reactor and accelerator data with the six oscillation parameters listed above 
is to identify the two mass splittings and the two mixing angles which drive the solar and 
atmospheric transitions with (Am| x , #12) and (|Am| 2 |, #23); respectively. The sign of the 
atmospheric mass splitting Am^ with respect to the solar doublet is one of the unknowns 
within the neutrino sector, i.e. we do not know if the neutrino mass spectrum is normal 
(Am 2 2 > 0) or inverted (Am| 2 < 0). The best fit point for the combined analysis of solar 
neutrino data together with KamLAND reactor data is at Am^ = 8.0 x 1CT 5 eV 2 
and sin 2 #12 = 0.31 3 . The 90% CL. allowed ranges of the atmospheric neutrino oscillation 
parameters obtained by the Super-Kamiokande experiment are |lz 



14. 



I Am|, I ~ lAmlil = (1.5 - 3.4) x 10~ 3 eV 2 , 0.36 < sin 2 23 < 0.64 (1) 

The mixing angle #13 (which connects the solar and atmospheric neutrino realms) and 
the amount of CP violation in the leptonic sector are undetermined. At present, the upper 

n 

bound on the angle #13 coming from CHOOZ reactor neutrino data [14J is: 

< sin 2 #13 < 0.04 (2) 

at the 90 % confidence level at Am^ = 2.5 x 10~ 3 eV 2 . This constraint depends on the 
precise value of Am^, with a stronger (weaker) constraint at higher (lower) allowed values of 
|Am§]J. Future reactor neutrino oscillation experiments could measure the value of sin 2 # 13 , 



1 We restrict ourselves to a three- family neutrino scenario analysis. The unconfirmed LSND signal cannot 
be explained in terms of neutrino oscillations within this scenario, but might require additional light sterile 
neutrinos or more exotic explanations ||J . The ongoing MiniBooNE experiment 

7] 

is expected to explore 

all of the LSND oscillation parameter space 



2 Am| = mf — Tfij throughoujLlhe paper 
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We use the notation of Ref. (llj throughout. 

ere, we will use |A 
Q and K2K Q 



4 For the numeric analysis presented here, we will use |Ato 2 9 | = 2.4 x 10 3 eV 2 , which lies within the best 



'321 

fit values for the Super-Kamiokande [12j and K2K [13( experiments. 
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as explored in detail in Ref. Q|. Current neutrino oscillation experiments do not have any 
sensitivity to the CP-phase 5. The experimental discovery of the existence of CP violation in 
the leptonic sector, together with the discovery of the Majorana neutrino character, would 
point to leptogenesis as the source for the baryon asymmetry of the universe, provided that 
accidental cancellations are not present. 

The main aim of this paper is a careful study of the sensitivity to the currently unknown 
parameters mentioned above, that is to the small mixing angle #13, to the ordering of the 
neutrino mass spectrum, and to the amount of CP violation in the leptonic sector, which 
could be achieved by exploiting the NuMI neutrino beamline. We thus concentrate on 
the NuMI beam potential exploited in an off-axis configuration as proposed by the NOz/A 



experiment 
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17| . The location of the far detector is at 10 km off-axis with a baseline of 810 
km. The mean neutrino energy is 2.3 GeV. We have considered two possible v e (z/ e ) detection 
techniques. First, the possibility of a 30 kton totally active low Z tracking calorimeter 
detector, as the one considered in the revised NOz/A proposal Q|. The efficiencies of such 
a detector for v e (p e ) identification is approximately 24% and the background is typically 
two-thirds from electron (anti)neutrinos in the beam produced from muon and kaon decays 
and one third from neutral current events faking electron neutrinos. An alternative detection 
method explored here is the one provided by a Liquid Argon TPC detector, as the technique 
described in the FLARE Letter of Intent 18]. The efficiency of such a detector to identify 
the u e (i> e ) CC interaction is 80% and the background is dominated by the intrinsic v e and 
v e components of the beam jisl ]. 

The statistics at the far detector is governed by the product of three parameters: the 
total number of protons on target (n p ), the detector mass (m D ), and the detector efficiencies 
to u e (u e ) identification (e). Here we have studied the physics potential of three different 
experimental scenarios: 

• Small: As a first step, we consider a "Small" experimental setup without a Proton 
Driver. The number of protons on target without a Proton Driver is 6.5 x 10 20 per 
year. We have considered three and a half years of running in each polarity (i.e. three 
and a half years of neutrino and three and a half years of antineutrino data taking). 
Consequently, in this first scenario, the statistical figure of merit n p x nip x e is equal 
to a total of 3.3 x 10 22 in units of number of protons times kton. This setup could 
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be achieved, for instance, with the 30 kton totally active low Z tracking calorimeter 
detector at 24% efficiency described above or with a 9 kton Liquid Argon TPC detector 
at 80% efficiency. 

• Medium: We study a possible upgrade of the Small configuration, referred to as 
the "Medium" experimental setup by increasing the statistical figure of merit by a 
factor of five. This statistics factor could be accomplished if the mass of the liquid 
argon detector is upgraded to 45 kton without a Proton Driver in three and a half 
years running in each polarity. Or, equivalently, the same statistics could be achieved 
in the NOvA. experiment running for four and a half years with a Proton Driver in each 
polarity (i.e. four and a half years of neutrino and four and a half years of antineutrino 
data taking). With a Proton Driver, the number of protons on target is assumed to 
be 25.2 x 10 20 per year. Thus the statistics running with a Proton Driver for four and 
a half years is five times the statistics for running without a Proton Driver for three 
and a half years 5 . 

• Large: The third experimental setup explored is a "Large" experiment, in which 
both the initial detector mass and the total number of protons on target are increased 
by a factor of five. The Large scenario could be obtained, for instance, by the 
combination of a 45 kt liquid argon detector with a Proton Driver running for four 
and a half years in both the neutrino and antineutrino modes. 

The ratio of the statistics in the three different experimental scenarios 
Small:Medium:Large is therefore 1:5:25. In the next section we review the neu- 
trino oscillation formalism, while numerical results are presented in the following sections. 



I. PRELIMINARIES 



Since we are exploiting the — > v e and z/ M — > v e appearance channels, the observables 
that we use in our numerical analysis are the number of expected electron neutrino and 
antineutrino events. For the central values of the already measured oscillation parameters, 

5 This factor of 5 is the ratio of the protons on target per year times the number of years of data taking 
with and without a proton driver: 2 6 5 5 2 ^ 3 4 5 5 . 
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we have thus computed the expected number of electron and positron events and 

at the far detector located 10 km off-axis at L = 810 km, assuming positive or negative 

hierarchies, which are given by: 

N^l% +) = J dE v L) a v(p) {E v ) P v ^ [9 ^ e) (E V} L, 13 , 5, 8 23 , 8 12 , ±Am^, Am 2 2l ) 

(3) 

where #23, #12 ; an d Am^ are taken as perfectly known, $„rp) denote the neutrino 

fluxes, and ov(z?) the cross sections. The neutrino (antineutrino) flux, which peaks at 2.3 GeV, 
is integrated over a narrow 1 GeV energy window (E min = 1.8 GeV and E max = 2.8 GeV). 

The Up — > u e ( —>■ u e ) appearance probabilities in long baseline neutrino oscillation 
experiments, assuming the normal mass hierarchy, read [3]: 

P„ M „ e = X + 8 2 + V + #cos(A 31 + 5) + P , 

P p , = X_8 2 - Y_8 cos(A 3 i - 5) + P . (4) 



In the last expressions, # = sin #13 and the coefficients X± and Y± are determined by 
I ZA 3 ism^a^ T ^3iJ I " 

^± = ^23 



(aL t A 3 i) 

v 4-0 nr~P~ +s J A 31 sin(aL^ A 31 )\ J A 2 isin(aL) \ 

V ± = ±2^P~ & = ±8c 12 s l2 c 23 s 23 I (qLtA3i) j ( QL j , (5) 



p _ r 2 .2 2fl / A 2 i sin (aL) " J 

^0 — C 23 Sm /y 12 S ^ 

where Ajj = \Am 2 j\L/4E, and a = GfN £ /\/2 denotes the index of refraction in matter, 
Gf being the Fermi constant and N e is a constant electron number density in the Earth. 
We denote the first, second and third terms in Eqs. (jlj) as the atmospheric, interference 
and solar terms, respectively. When #13 is relatively large, the probability is dominated 
by the atmospheric term. Conversely, when # 13 is very small, the solar term dominates. 
The interference term is the only one which contains the CP phase 5, and it is clear from 
Eqs. flU) that it is also the only one which differs for neutrinos and antineutrinos besides 
matter effects. 

We can ask ourselves whether it is possible to unambiguously determine #13 and S by 
measuring the transition probabilities z/ M — > u e and —>■ v e at fixed neutrino energy and at 
just one baseline. The answer is no. At fixed neutrino energy E v and baseline L, if (# 13 , 5) 
are the values chosen by Nature, the conditions: 



can be generically satisfied by another set {0' 13 , 5'), known as the intrinsic degeneracy [19] . 
It has also been pointed out that other fake solutions might appear from unresolved degen- 
eracies in two other oscillation parameters: 

1. The sign of the atmospheric mass difference Am^ may remain unknown. In this 



particular case, P(0' 13 , 5', — Am| 2 } 



P(9 13 ,S, Ami 



32 I 



22 



231 ] . More specifically, the 



v e ) appearance probabilities for the inverted hierarchy are 



P 



X_9 2 + Y_6 cos(A 31 - 5) + P e 



X+9 2 - Y+9 cos(A 3 i + 5) + P Q . 



(6) 



2. Disappearance experiments only give us information on sin 2 2^23: is #23 in the first 
octant, < #23 < 7r/4, or is it in the second one, 6*23 — > 7r/2 — 6*23?. In terms of 



oscillation probabilities, P( 



^13' ^ ) 2 



023) = p(0i3,s,e 23 ) mm 



Extensive work has been devoted recently to eliminate such fake solutions. In simple 
terms, when considering data from two or more experiments, degenerate solutions may occur 
at different locations in parameter space for different experiments, and therefore could be 



excluded 
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25l . l2q . As shown in Ref . 



there exists a simple way to understand if 



the various degeneracies are eliminated when several experiments are combined. We briefly 
review here the analysis of Ref. Q], and illustrate its results for the NuMI beamline exploited 



in an off-axis mode. In the overlap region of the — > v e bi-probability diagram |2lJ, see 
Fig. CQ)(a), there exist generically four solutions 6 for the unknown oscillation parameters 
#13 and 5. Two solutions correspond to the normal hierarchy Q| and have approximately 
equal values of sin 5 but different values of the sign of cos 5. The other two solutions are for 



the inverted hierarchy 



3, 



231 ] . and they also have approximately equal values of sin 5. This 



value of sin S for the inverted hierarchy is different than the value of sin 5 for the normal 
hierarchical case. In Ref. (2^, an identity connecting the difference between the mean 



6 We assume 6*23 = 7r/4; we will show that small variations of this parameter do not affect the results 
presented here in a significant way. 
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NuMI: E = 2.3GeV and L=810km 



1.0 



NuMI: 10km <E>=2.3GeV L=810km 
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FIG. 1: (From ref. j^)/ ) (a) The bi-probability diagram for NuMI 10 km off-axis showing the allowed 
regions for both the normal (dashed) and inverted (dot-dashed) hierarchies as well as the ellipses 
/or sin 2 2^i3 = 0.05. The large "-)-" marks the neutrino and anti-neutrino probabilities with the CP 
phase, 5 = vr/4, assuming the normal hierarchy. The ellipses and point along the diagonal labeled 
critical correspond to the largest values for which there is overlap between the normal and inverted 
hierarchies. (b)The allowed regions in the sin 2 2#i3 v sin 5 plane for the NuMI 10 km off-axis 
experiment, assuming the true solution is the normal hierarchy with sin 2 2#i3 = 0.05 and 5 = 7r/4 
("-\-"). The upper blue (lower red) contours are for the normal (inverted) hierarchy whereas the 
solid (dashed) contours are for cos 5 > (< 0). The experimental setup corresponds to the Large 
scenario considered here. The ellipses correspond to 68, 90 and 99% C.L. contours. 

values of sin 5 (which governs the amount of leptonic CP violation) for the two hierarchies, 
to the mean values of #13 for both hierarchies, is derived. Such an identity turns out to be 
extremely helpful in understanding if the combination of several experiments can eliminate 
the fake solutions, since the location of the fake solutions in the (sin 2 2^13, sin S) plane can 
be computed in a straightforward manner. If we apply this identity relating the solutions 
corresponding to the positive and negative hierarchy to the NuMI 10 km off- axis experiment, 
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it was found: 

<sin5> + -(sin5>_ = 1.41^^, (7) 

where (sin £)+(-) are the mean values of the two solutions of sin 5 for each hierarchy, see 
Ref. |20j for details 7 . For the sake of illustration we show in Fig. (fl|)(b) the \ 2 contours 
for NuMI 10 km off-axis in the Large experimental setup explored in the present study, 
assuming that the true solution is the normal hierarchy and that the values of (sin 2 2#i3, 
5) are (0.05, vr/4), respectively. The NuMI 10 km off-axis is operated above oscillation 
maximum: there are thus four solutions in the (sin 2 2# 13 , sin S) plane. 

The existence of such a simple relation, Eq. ((7j), among the true and fake solutions in 
terms of sin 5, together with the fact that it is precisely sin S the quantity which drives 
the amount of leptonic CP- violation, has motivated us to consider sin 5 as the relevant 
parameter in our analytical and numerical studies. In the case of the T2K experiment, the 
difference between the true and fake solutions for the CP violating parameter sin 5 is 0.47 
at sin 2 2#i3 = 0.05. This factor of 3 decrease with respect to the NuMI off-axis experiment 
is primarily due to the T2K baseline is 1/3 the NuMI off-axis baseline. 

If the mixing angle 8 2 s ^ vr/4, the (sin 2 26 , 13 , sin S) plane should be translated into the 
(2 sin 2 9 2 3 sin 2 26>i 3 , y/2 cos #23 sin 5) plane, as shown in Ref. (01) ■ Assuming this mapping, 
the results presented in the next sections would be almost identical even if 6*23 7^ 7r/4. 

II. sin 2 26»i3 SENSITIVITY 

In the present section we explore the sensitivity of the proposed NuMI long baseline 
off-axis experiment to a — + v e {v^ — > v e ) oscillation search in the appearance mode. In 
Figs. (J2J), © and (jlj we depict the sensitivity contours in the (sin 2 2#i3, sin S) plane for the 
three different setups described in the introduction for neutrino data, antineutrino data and 
for a combined analysis of both neutrino and antineutrino data, respectively. 

We have four different curves for neutrinos (for antineutrinos as well as for the combi- 
nation of the two channels), according to the mass spectrum hierarchy and to the sign of 

7 If the detector is located 12 km off-axis at the far site, the numerical factor in Eq. Q) is 1.46 instead of 
1.41. Consequently, the changes associated with placing the detector at 12 km off-axis rather than at 10 
km off- axis are small. 



s 



cos S. Since both the sign of the atmospheric mass splitting and the sign of the cos 5 may 
remain unknown, the maximum sensitivity to sin 2 2#i3 versus sin 5 in a given setup should 
be identified with the most conservative curve among the four possibilities. We describe in 
detail Fig. (0) but the same criterion should be applied to Figs. (j2J) and (@J). In the left panel 
of Fig. (J2J) it is depicted the sensitivity to sin 2 2# 13 versus the CP violating quantity sin 6 
in the two conservative pictures, i.e. when the sign of cos 5 < in the normal hierarchical 
scenario or when cos 5 > in the inverted hierarchy picture. Being the sign of cos S and 
the neutrino mass hierarchy unknowns within the neutrino mixing sector, the sensitivity 
to sin 2 6*i3 must be associated with the tightest bound among the two possibilities. In the 
particular case that we are describing here, i.e, the one exploiting the neutrino channel in- 
formation, the sensitivity curve is given by the red solid curve in the left panel of Fig. (j2J) . 
In the right panel of Fig. (J2J) we show the most optimistic scenarios where cos 5 > in the 
normal hierarchy picture or cos 5 < in the inverted one. However, we should remark here 
that these curves do not represent the true sensitivity to sin 2 26 13, which is the first priority 
of the near future neutrino oscillation experiments, before the measurements of the sign of 
the atmospheric mass difference and the sign of cos S. 

We point out here the existence of zero-mimicking solutions: there exists, at a fixed 
neutrino energy and baseline, a point in the (sin 2 2^13, sin 5) plane at which the first (at- 
mospheric) and second (interference) terms in Eq. (j3J) exactly cancel, and the situation is 
indistinguishable from the one in which sin #13 = 0. In vacuum, the location of the zero- 
mimicking solution at a given S is 

sin 26 13 = _ 2 _^^li cos(A32 ± s) f (8) 
smA 31 tan t) 2 3 

where the sign +(-) refers to neutrinos (antineutrinos). If the experiment were operating at 
the vacuum oscillation maximum, the zero-mimicking solution would be located at: 

sin 2#i2 „ , x 

sin 26 13 = ±2 A 21 — — - sin 5 , (9) 
tan c?23 

where the sign +(-) refers to neutrinos (antineutrinos). At the vacuum oscillation maximum 
the sin S value of the zero-mimicking solution is positive for neutrinos (negative for antineu- 
trinos). Figure ©(b) depicts the zero- mimicking solution for the T2K experiment (27 1. 
T2K will use a steerable neutrino beam from JHF to Super-Kamiokande and/or Hyper- 
Kamiokande as the far detector (s). The mean energy of the neutrino beam will be tuned 



sign(Am| 1 )cos 6 < sign(Am| 1 )cos 6 > 




sin 2 2(9 13 sin 2 2(9 13 

FIG. 2: T/ie 99% CL sensitivity contours in the (sin 2 2#i3 ; sin<5,) plane for sin 2 26 '13, exploiting only 
the data in the neutrino channel. In the left panel we depict the sensitivities for the most restrictive 
choice of parameters, which corresponds to cos 5 < and normal hierarchy ( dashed blue curve ) or 
to cos 5 > and negative hierarchy (solid red curve). In the right panel we show the sensitivity for 
the choice of parameters where cos 5 > in the normal hierarchy (solid blue curve) or cos 5 < 
and the hierarchy is inverted (dashed red curve). The labels L, M and S correspond to the Large, 
Medium and Small experimental setups explored in this study, respectively. See Table in the 
Appendix for numerical limits. 

to be at vacuum oscillation maximum, A 31 = |, which implies a mean neutrino energy 
(E u ) = 0.6 GeV at the baseline of 295 km, using lAmf-J = 2.4 x 10~ 3 eV 2 . This neutrino 
energy can be obtained with a 3° off-axis beam. Since T2K will be operated at vacuum 
oscillation maximum, this experiment is insensitive to the CP conserving quantity cos<5. 

When matter effects are considered, the situation is slightly more complicated. The zero- 
mimicking solution is in general located at different points in the (sin 2 26 l i3, sin 5) plane, as 
illustrated in Fig. ©(a), showing the zero-mimicking solution for the NuMI 10 km off-axis 
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sign(Am| 1 )cos 6 < sign(Am| 1 )cos 6 > 




sin 2 2(9 13 sin 2 2(9 13 



FIG. 3: The 99% CL sensitivity contours in the (sin 2 2#i 3; sin <5 ^ plane for sin 2 28 'i 3 , exploiting only 
the data in the antineutrino channel. The curves and labels are same as those of Fig. 

experiment 8 : 

• oa sin2g 12 ( Agi sin (aL) | / (aL =f A31) | . f . 

sin 20i3 = -2- — — 4 = f <^ -r . . r ^ cos A 32 ± 5) , 10 

tan0 2 3 I J I A 3i sm(aL =f A 31 ) J 

The existence of zero-mimicking solutions allow us to understand the shape of the sensi- 
tivity curves, i.e. the Figs. (j2J), (JHJ) and (jlj). Exploiting the neutrino data, the sensitivities 
will therefore improve enormously as the experimental setup (i.e. the statistics) is improved, 
as long as sin 5 is negative. The sensitivity to sin 2 2# 13 is maximal when sin 5 = —1. On 
the other hand, if the antineutrino data is exploited, the situation is reversed: the sensitiv- 
ities will improve in a significant way as the setup is upgraded if sin 5 is positive, and the 
optimal sensitivity in this case is reached when sin 5 = 1. When combining the data from 
the neutrino and antineutrino channels, the sensitivity curves are flatter than in each sepa- 
rate case (i.e. considering only the neutrino data or only the antineutrino one), see Fig. 
This flattening effect on the sensitivity curves when adding the information of both channels 



The NuMI beam energy (2.3 GeV) is about 30% above the vacuum oscillation maximum energy for its 
baseline, i.e 810 km. Since NuMI 10 km off-axis is operated above oscillation maximum, this experiment 
is sensitive to the sign of cos 5. 
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sign(Am| 1 )cos 6 < sign(Am| 1 )cos 6 > 




sin 2 2(9 13 sin 2 2(9 13 



FIG. 4: T/ie 99% sensitivity contours in the (sm 2 2#i3 ; sinS) plane for sin 2 28 '13, exploiting both 
the data in the neutrino and in the antineutrino channel. The curves and labels are same as those 
of Fig. 0J. 

increases as the exposure decreases: in the Small experimental setup considered here the 
curves are basically flat in the (sin 2 26 l i3, sin 6) plane, see Fig. (j3J). 



III. MASS HIERARCHY DETERMINATION 



In this section we study the possible extraction of the sign of the atmospheric mass 
splitting with the NuMI 10 km off-axis experiment, within the context of the three reference 
experimental setups considered in the present paper. We have performed a x 2 analysis of 
the data in the (sin 2 20 13 , sin 5) plane. We assume Nature has chosen the normal or inverted 
hierarchy and we attempt to fit the data to the expected number of events for the opposite 
hierarchy. Generically one expects two fake solutions associated with the wrong choice of 
the hierarchy at fixed neutrino energy and baseline. The x 2 function, from the combination 
of the neutrino and antineutrino channels, reads: 

^E(^p) 2 . (") 



p=e^,e~ 
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FIG. 5: (a) Zero-mimicking solution for the NuMI off-axis experiment, in the fsin 2 2#i3, sin<5j 
plane, for neutrinos (upper curves) and antineutrinos (lower curves). The blue (red) depicts normal 
(inverted) hierarchy. The vacuum zero-mimicking solutions (not depicted here) are located between 
the two hierarchies, therefore matter effects are small. The solid (dashed) curves depict the positive 
(negative) sign for cos 5. (b) The same as in (a) but for the T2K experiment, as an illustration of 
the zero-mimicking solutions at vacuum oscillation maximum, A31 = 5. 

where 5J\fp is the statistical error on Mp , the simulated data 

JV+ = Smear (N+ + B p ) , (12) 

where B p is the number of background events when running in a fixed polarity p and we have 
performed a Poisson smearing to mimic the statistical uncertainty. In Fig. © we depict the 
results for the sign(Am3 1 )-extraction by exploiting the neutrino and antineutrino data in the 
three reference setups. As expected, the best sensitivity is reached with the most ambitious 
scenario, relying on the proton driver option. The shape of the exclusion lines can be easily 
understood in terms of matter effects, which are quite significant for the NuMI off-axis 
experiment as can be clearly noticed from the bi-probability diagram, Fig. ©(a). The shift 
observed in the bi-probability events is proportional to the size of the matter effects, which 
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0.02 0.04 0.06 0.08 0.1 
sin 2 20 13 



FIG. 6: Sensitivity to the sign (Am^j- extraction at the 95% CL within the three reference setups 
explored in the present study. The labels L, M and S correspond to the Large, Medium and 
Small experimental setups explored in this study, respectively. The dashed black curve is obtained 
from Eq. setting (sin5)_ = —1 (^(sin5) + = +1) for the normal (inverted) hierarchy. This is the 
bound that would be obtained with infinite statistics and in the absence of backgrounds. 



are obviously crucial to resolve the hierarchy of the neutrino mass spectrum 9 . The sensitivity 
to the measurement of the sign of the atmospheric mass difference is expected to be better 
when the sign of sin<5 is negative: in the case of the Medium experimental setup, the 
sensitivity to the sign (Am^-extraction is lost for positive values of sin 5. We show as well 
in Fig. (jUJ) the theoretical limit on the sign(Am| 1 )-extraction, which acts as a rigorous upper 
bound on the experimental sensitivity curves. A possible way to resolve the fake solutions 
associated to the sign of the atmospheric mass difference would be to combine the data from 
the proposed NuMI 10 km off-axis and T2K experiments [^J 2^j]. The complementarity of 
the NuMI and T2K experiments can be explicitly shown by exploiting the identity given in 



9 Recently, new approaches for determining the type of hierarchy have been proposed 28] by exploiting other 
neutrino oscillations channels, such as muon neutrino disappearance, and require very precise neutrino 
oscillation measurements. 
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the introductory Section by Eq. (JJJ) and in Ref . 20] . The difference in the location of the 
fake solutions associated to the wrong hierarchy in the (sin 2 26 l i3, sin 5) plane for these two 
experiments reads: 



(sin 5) 



T2K 
fake 



(sin 5) 



NuMI 
fake 



0.94 



sin 2 29 
0.05 



13 



(13) 



This relation implies that the wrong solutions would appear in different regions of the 
parameter space for the two experiments, and therefore the fake solutions could be eliminated 
by a combined NuMI/T2K analysis. 10 

IV. CP-sin 5 MEASUREMENT 

In the present section we explore the sensitivity to CP violation for the three different 
experimental scenarios under consideration. The results are summarized in Fig. (J7J), in 
which we depict the exclusion contours at the 95% CL corresponding to the Large and 
to the Medium setups. The measurement of leptonic CP violation is certainly not within 
reach for the less ambitious scenario, i.e. the Small experimental setup described in the 
introductory section, not shown in Fig. (J7|). 

The exclusion lines depict the value of sin 5 at which the 95% CL error in the CP violating 
parameter sin 5 reaches sin 5 = 0, corresponding to the CP conserving case. An important 
point to note here is that, in order to compute the curves in Fig. (JJJ), we have considered the 
impact of the degeneracies associated to the sign of the atmospheric mass difference. It may 
happen that the error on the fake-sign solutions reaches sin 5 = and therefore the analysis 
would be consistent with CP conservation even if the error in the true solution has not 
reached sin 5 = 0. The mass hierarchy-sign degeneracies affect the CP-sensitivity contours 
only if sign(Am3 1 ) sin 5 is positive, as would be expected from the analysis performed in the 
previous section to the sensitivity to the sign of the atmospheric mass difference. 



Based on previous work [25|, the authors of Ref. 29] have recently shown that it would be possible to 
extract the sign of the atmospheric mass difference by exploiting the NuMI off-axis beamline with just a 
neutrino run, provided that two detectors would be placed at the same E/L. The experimental picture 
would be given by a near detector located before the NOtA far site (probably at 200 km, to optimize the 
sensitivity) and a second detector at the far site. 
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FIG. 7: Sensitivity to CP violation (i.e. sin 5 ^ 0) at the 95% CL. The difference in the shape 
of the two top exclusion curves is due to the additional restrictions coming from the presence 
of the signfAm^J-degeneracies. If the hierarchy were determined, the top curves would be close 
to a reflection of the bottom ones. The labels L and M correspond to the Large and Medium 
experimental setups explored in this study, respectively. 

V. ANALYSIS IN ENERGY BINS 

In the present section we restrict the analysis to the intrinsic degeneracies: we assume 
that Nature has chosen Am^ > and #23 — 7r /4, and we concentrate on only one scenario, 
the Large setup. 

In order to eliminate the fake solutions associated with the wrong choice of the sign of 
cos S, we exploit here the energy dependence of the signal. We have thus divided the total 
number of events in two bins of equal width AE U = 0.5 GeV 11 . The \ 2 function of Eq. 
at the fixed baseline L = 810 km now reads: 



A very conservative estimate for the neutrino energy resolution in the range of interest is AE V / E v ~ 50% 




(14) 
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Assuming a positive sign for the atmospheric mass splitting, the result of the fit assuming 
no energy binning in the signal for a particular central value in the (sin 2 2#i3, sin 5) plane is 
shown in Fig. (J8J)(a). where two solutions arise, one of them (the fake one) corresponding 
to the wrong sign of cos 5. If the data analysis is performed with the opposite sign of the 
atmospheric mass splitting (negative), the two fake solutions associated to the wrong choice 
of the spectrum hierarchy are not present at the 95% CL for the particular central value 
chosen in Fig. (JHJ). 

If one exploits the energy information in the signal by performing an analysis in energy 
bins, the intrinsic, fake solution is resolved, as it is shown in Fig. (jEJ)(b). It has been shown 
that, for sufficiently large 6*13 and in the vacuum approximation, apart of the true solution, 
there is a fake one at 19j: 

sin 5 ~ sin 5 , 

6> 13 ~ #13 + cos 5 sin 2^12 A 2 i cot #23 cot (A31) . (15) 

The fake solution would be located at a value of sin 2 29' 13 which is energy-dependent. The 
degeneracies associated with the different energy bins would therefore have different loca- 
tions in sin 2 26^3. We illustrate the results for the two bins separately in Fig. (|9*|). The 
fake solutions appear in two different regions of the parameter space: when combining the 
information from the two separate bins, the intrinsic degeneracy would disappear. In order 
for this conclusion to hold true, the energy dependence of the signal has to be significant 
enough; otherwise, the analysis in energy bins would not provide an effective elimination 
of the fake solutions. We have performed the analysis in energy bins for smaller values of 
sin 2 2# 13 than the one shown in Figs. (jSJ) and (|9*|). We find that the intrinsic degeneracy is 
resolved if sin 2 2#i 3 > 0.02. 



VI. CONCLUSIONS 



We show the neutrino oscillation physics potential that can be achieved with the Fermilab 
NuMI beamline 10 km off-axis and three different experimental setups, differing in the proton 
luminosities and/or in the detector sizes. We provide a complete study of the sensitivities to 
sin 2 2#i3, to the hierarchy of the neutrino mass spectrum, and to the CP violating parameter 
sin S for the three different scenarios. We present our results in the (sin 2 28is, sin 5) plane; 
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0.02 0.04 0.06 0.02 0.04 0.06 

sin 2 20 13 sin 2 20 13 

(a) (b) 

FIG. 8: (a) 95% CL contours from a simultaneous x 2 fit to sin 2 2#i3 and sinS without binning the 
data. The central value is denoted by a star: sin 2 2#i3 = 0.03 and sin 5 = —0.5. We have included in 
the former analysis statistical errors and backgrounds. The true (fake) solutions associated with the 
CP violating parameter sin 5 are depicted in solid (dotted) contours. (b)The result of the combined 
X 2 analysis of the two energy bins. 

this choice helps enormously in understanding the location of the solutions for different 
experiments, and turns out to be very easy to generalize if 8 2 3 7^ 7r/4. We also explore 
the benefits of a modest energy resolution: with a 50% energy resolution, the intrinsic 
degeneracies are lifted if sin 2 2# 13 > 0.02. 
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FIG. 9: 95% CL contours from a simultaneous \ 2 fit to sin 2 2#i3 and sin5 applying a binning to the 
data. The results of the analysis with the first (second) bin of data are depicted in blue (magenta). 
We have included in the former analysis statistical errors and backgrounds. The fake solutions 
disappear when the two bins of data are combined, as shown in Fig. f$)(b). 



Energy. 
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VIII. APPENDIX 



A. Oscillated statistics 

In Table (JJ) we provide the computed charged-current event rates at the NOi/A far site 
(810 km) in the Medium experimental setup described in the Introductory Section. 



v e (signal) 


v e (background) 


v e (background) 


Ve (013 = 0) 


145 


50.0 


2.87 


7.55 


v e (signal) 


v e (background) 


9 e (background) 


Ve (013 = 0) 


44.8 


6.64 


17.4 


2.33 



TABLE I: Calculated charged currents neutrino and antineutrino event rates (signal and back- 
grounds) for NuMI (baseline of 810 km, 10 km off-axis) in the Medium experimental setup de- 
scribed in the Introduction. In order to compute the signal we have assumed energy independent 
oscillation probabilities P (v^ — ► v e ) and P (v^ — ► v e ) equal to 1%. The v e (Qy& = 0) and v e 
(#13 = 0) are the contributions from P & in Eq. The non oscillated (v^) event rate can be 
computed by multiplying the v e ( D e ) signal by 100 and dividing the result by the v e ( v e ) detection 
efficiency. 
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B. sin 2 2#i3 sensitivity 





mode/Setup 


Small 


Medium 


Large 


sin 2 20i3 
99% CL-sensitivity 
(most restrictive) 


V 


2.0 x 1(T 2 


1.6 x 10~ 2 


1.0 x 10~ 2 


V 


4.0 x 10~ 2 


2.0 x 10" 2 


1.3 x 10~ 2 


V + V 


1.2 x 10~ 2 


7.0 x 10" 3 


4.0 x 10~ 3 


sin 2 20i3 
99% CL-discovery 
(least restrictive) 


V 


4.0 x 10" 3 


1.2 x 10" 3 


3.0 x 10" 4 


V 


1.0 x 10~ 2 


4.0 x 10~ 3 


1.5 x 10~ 3 


V + V 


4.0 x 10~ 3 


1.2 x 10~ 3 


3.0 x 10~ 4 



TABLE II: 99% CL sensitivity limits to sin 2 20i3, extracted from Figs. HJ) ; |3J) and using only 
the neutrino channel, only the antineutrino one, or both channels, respectively. The first three 
rows indicate the most restrictive limits which is independent of the value of the CP phase 5 and 
of the type of hierarchy chosen by nature. The last three rows show the 99% CL-sin 2 20i3 discovery 
limits for the most favorable choice of the parameters. The three columns refer to the three different 
experimental setups explored here. 
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